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Abstract: We consider steady states for a class of mechanical systems with particle- 
disk interactions coupled to two, possibly unequal, heat baths. We show that any steady 
state that satisfies some natural assumptions is ergodic and absolutely continuous with 
O respect to a Lebesgue-type reference measure and conclude that there exists at most one 

^L, absolutely continuous steady state. 

^ . 

Introduction. 

Explaining macroscopic phenomena of open systems from the microscopic dynamics 
is an intriguing subject in statistical mechanics. By an open system we mean a deter- 
ministic and energy conserving system that exchanges energy and matter with multiple 
heat baths. Some examples of such systems have been studied |[T1l2l[^ffll5ll6ll7l[8l l9lfT0l 
fm . but a good understanding of non-equilibrium behavior, i.e. when the heat bath are 
not equal, has not been reached yet. 

In this paper we consider a mechanical system that attempts to reproduce the phe- 
nomenological laws of thermodynamic transport [9 |. It consists of a rectangular domain 
containing a chain of N identical pinned down disk scatterers that are allowed to rotate 
freely. See Figure Q] The particles in the system bounce elastically from the walls and 
exchange energy with the disks through "perfectly rough" collisions ||8]|9). The system 
is coupled to two possibly unequal heat baths through the openings, which correspond 
to the "left" and the "right" sides of the rectangle. Once a particle reaches an opening, 
it leaves the system forever. New particles can be introduced to the system by the heat 
baths; they are emitted at random times according to some probability distributions that 
describe injection positions and velocities. Neither energy nor the number of particles 
is conserved, so the system can be loosely referred to as a "grand-canonical ensemble". 
The detailed settings are available in section Q] 

This paper demonstrates rigorous results regarding invari ant densities and ergodic- 
ity of the steady states, which we refer to as invariant measures. Theorem Q] provides 
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Fig. 1. ID array of disks in a rectangle 



conditions under which an invariant measure is guaranteed to be ergodic and absolutely 
continuous with respect to a natural reference measure. From Theorem [T] we immedi- 
ately obtain uniqueness of the absolutely continuous measure as well as the ergodic 
decomposition. In the equilibrium situation, an absolutely continuous invariant mea- 
sure can be written down explicitly [7|. By Theorem[T]we conclude that this invariant 
measure is ergodic. 

The question of existence and uniqueness of the steady states for other types of 
open systems that attempt to model thermodynamic transport has been studied to some 
extent. The best studied examples are chains of anharmonic oscillators coupled at both 
ends to heat reservoirs [11]. The existence and uniqueness results for such systems 
have been obtained in [ 3 5 , 6 1 under rather restrictive assumptions on the potential in 
the chain and coupling to the reservoirs. 

For open particle systems with indirect particle interactions, existence of the sta- 
tionary states has only been demonstrated in equilibrium situations by providing them 
explicitly [1,7]. For non-equilibrium states, existence turns out to be a nontrivial mathe- 
matical problem which requires a detailed understanding of the dynamics and we leave 
it for future work. 

Uniqueness results for open particle systems are scarce. Ergodicity and thus unique- 
ness of natural invariant measures for a 1-dimensional particle model was obtained in 
(TJ; we borrow some of the ideas developed there in our proof. For planar geometries, 
an important first step towards showing ergodicity we rely on was done in H, where 
the authors demonstrated that the action of the baths can drive the system from "almost" 
any state to "almost" any other state in a finite time. We chose simpler geometry in our 
model in order to focus on the essential properties of the system that lead to absolute 
continuity and ergodicity of the invariant measures leaving geometric complications 
aside. However, combined with the results in [4|, only minor modifications are needed 
for our argument to carry through for the class of systems described there. Our argu- 
ment might potentially be generalized to a wider class of geometries, e.g. presented in 

GEO- 

The organization of our paper is the following. Section Q] contains a detailed de- 
scription of the model. We state the main results in section [2] The remaining part of 
the paper is devoted to the proof of Theorem [1] We start with an outline of proof in 
section[3j where we state three propositions and show how they imply TheoremQ] The 
propositions are proven in sections |4] E] |6] andQ 



1. Model Description. 

1.1. Dynamics of closed systems. Let To be a rectangle bounded by y = ±1, x = 
and x — 2N, where N is an arbitrary positive integer. In the interior of Jq lie N disks 
Dj, 1 < j < N, of equal radii R < 1 centered at (2j -1,0) G T , 1 < j < N. 
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The centers of the disks are fixed and the disks are allowed to rotate freely around the 
center, each carrying a finite amount of kinetic energy derived from its angular velocity. 
Denote the states of the disks by ((pj,ojj), 1 < j < N, where ipj is the disk's angular 
position relative to a marked reference point and uij is the disk's angular velocity. 

A number of particles move around in the playground r = To \ yj^ =l Dj, with par- 
ticle positions qi E T and velocities Vi E M 2 . Apart from collisions with the boundary 
of the playground dT = <9-To U (UjL 1 dDj) the particles are assumed to move freely 
with constant velocities; particles do not interact with each other. The collisions with 
dTo are specular, and upon collision of a particle with a disk, a certain energy exchange 
occurs |89|. More precisely: 

The phase space of such system with k particles is 

fi k = (r k x dDi x • • • x dD N x M. 2k+N )/ ~, 

where 

q = {q\, ■ ■ ■ ,qk)Er k denotes the positions of k particles, 

cp = (ipi, • • • , <pjy) E dDi x • • • x 8Dm denotes the angular positions of the disks, 

\r = (vi, ■ ■ ■ ,Vk) E M. 2k denotes the velocities of k particles, 

ui = (u)i , • • • , ojjv) E M. N denotes the angular velocities of the N disks, 

and ~ is the relation identifying pairs of points on the collision manifold: Mf. = 

{(q, ip, v, ui) : qi E dr for some i} with the rules of identification as follows: 

Let Vi — ((vi)t, (vi)±) be the tangential and the normal components of Vi. If a par- 
ticle collides with the boundary of the playground, qi E Zo, then the angle of reflection 
is equal to the angle of incidence, i.e. 

= -(v%)±, (Vi)'t = (Vi)t 
If collision with a disk occurs, E dDj for some j, then 

= ~(vi)x, (Vi)' t = (Vi) t - ——((vi)t - R0Jj), 

I + 77 

2 

Rw'j = Ruj + - — {{vi) t - Ruj), 

where 77 = ^rj- is a dimensionless parameter relating the moment of inertia of the disc 
(9, the mass of the particle m, and the radius of the disc R [8 ,9|. Throughout the paper 
we assume that < 77 < 00. 

The identification ~ in the definition of is as follows: X ~ X', X, X' E Mk, 
if the coordinates of X and X' are equal except for Vi's such that qi E dr. If qi E Zo, 
we replace v.- L in X by v[ in X' . If ^ E dDj, we replace Vi and uij in X by v[ and 
ujj in X'. Note that simultaneous collisions of several particles with the same disk are 
not defined, while there is no problem with simultaneous collisions with different disks 
and/or OTq. 

Define the discontinuous flow <P T on ilk by 

<P T (q, ip, v, ui) = (q + vr, <p + lot, v, ui) 

if no collisions are involved. When collisions occur, the rules of identification are given 
above. Then mk is an invariant measure for <P T on where 

m k = (A 2 |r)' £ x p\ dDl x • • • x p\ d D N x X 2 k+N- 

Here is d-dimensional Lebesgue measure and p\dD is the uniform measure on the 
circle dDj Q. 
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1.2. Dynamics of open systems. 

1.2.1. Coupling to heat baths. Suppose the rectangle Zo has two openings, = 
{0} x [-1, 1] £ dr a and 7^ = {2iV} x [-1, 1] e <9r , each connected to a heat 
bath that absorbs and emits particles. A particle absorbed by one of the baths leaves the 
system forever. The injection process for each bath is characterized by the following 
parameters: 

ql and qr - the injection rates of the baths. The injection processes are Poisson with 
rates ql and qr respectively. 

Tl(v) and Tr{v) - the distributions of the positions of injection with values in and 
7r- 

Al(Q and Z\l(C) - the distributions of the angles of injection with values in (— ~). 
Sl (0 and <Sr (0 - the distributions of the injected particle speeds with values in (0, 00). 

i.e. a particle is injected from the left bath at a random time r e (0, 00) given by 
exponential distribution of rate ql, with random position ( e 71 drawn from Yl(v), 
at random angle 5 drawn from Al(Q and with random speed s drawn from Sl($). 
Similarly for the right bath. 

We assume that each of the distributions Yl(v), Yr(v), Al{Q), Ar(£), Sl{$), and 
Sr(s) has positive density on the specified domains. 

7.2.2. The phase space. Since we are interested in the invariance properties, we would 
like to treat all particles in the system as identical and indistinguishable. Let Qk be as in 
subsection ll.il Then the phase space of the system coupled to heat bath(s) is a disjoint 
union 

to = ^T=otok, 

where [2% is the quotient of obtained by identifying the permutations of k particles. 
If we denote unordered sets by {...}, under phase space, points in J? are denoted by 

X = ({<?iv ,qk},(<Pi,-" ,Vn),{ v i>-" ,Vk}, (wi, • • • ,ujv)), 

with Vj understood to be attached to qj. Denote the quotient of the measure rhk by mk- 
Let <P T be a continuous-time Markov process on Q defined as follows: 

- <P T is the quotient of <P T identifying the permutations while no particles enter or exit 
the system; 

- if <P T € tok for all < r < tq, and a particle exits the system at time tq, then <P To [z) 
jumps to J?fe_i; 

- if ^ T € fik for all < r < tq, and a particle is injected from one of the baths or 
point sources at time tq, then <P To jumps to Qk+i- 

It is hard to write down the transition probabilities of <P T explicitly because particles 
may enter at any time. Denote the time-r transition probability starting at state X for 
<P T by if defined. PJr if not defined for all r and X: if, for example, X is a state 
such that two particles will have their first collision at the same time t with the same 
disk Dj, then PJ^ is not defined for all r > t. Note that is defined if the dynamics 
is defined with probability 1 . 
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Definition 1. We will say that a measure /i is absolutely continuous with respect a mea- 
sure v if for any measurable set A, v(A) = => n(A) = 0. We will say that two mea- 
sures [i and v defined on a set C are singular if there exist sets A and B, AC\ B = 0, 
A U B = C such that v(A) = and fi(B) = 0. 

When a measure /i is absolutely continuous with respect to v, we will denote it by 
fi <C v; we will also say that [i has a density with respect to v. If /i and v are singular, 
we will denote it [i _L v. 

Definition 2. Let m be a measure on Q such that for each nonnegative integer k, the 
conditional measure on [2^ is tti^. 

For our purposes, to is a natural reference measure: if p is any measure on Q, we 
are interested whether /i is absolutely continuous with respect to m. In the rest of this 
paper, we will refer to m as Lebesgue measure and when we will say that p is abso- 
lutely continuous without mentioning any reference measure, we would mean that [i is 
absolutely continuous with respect to m. 

1.3. Problems of interest. Given a measure v, the push forward of v under <P T , if de- 
fined, is given by (4> T )^v = j Q {PJ c ){A)dv{X). Note that P£ = (0 T )*<5 X when de- 
fined. 

Definition 3. A Borel probability measure [i on Q is called invariant under <P T if its 
push forward under <P T is defined for all r and 

fjt(A) = (($ T ).H)(A) - / (P T x )(A)dp(X). 

A classical way of analyzing flow <P T is studying the properties of its invariant mea- 
sures such as existence, uniqueness, absolute continuity with respect to the Lebesgue 
measure to, and ergodicity. Proving the existence of invariant measures would require 
very technical arguments that involve dealing with tightness and discontinuities; we 
leave these problems for future work. In this paper we are going to focus on the is- 
sues of absolute continuity with respect to to and ergodicity of the invariant measures 
provided they exist. 

The main result of this paper claims that, for the class of systems in consideration, 
if there exists an invariant measure p such that the measure of the set of all states with 
"trapped" particles is zero, then [i is both absolutely continuous and ergodic. 

2. Resuits. 

Definition 4. A state X e fl is said to contain a trapped particle if either 

— the velocity of the particle is zero, v — 0, or 

- the x-component of the velocity of the particle is zero, v x — 0, and the position of 
the particle has its x-coordinate "between the disks", i.e. 

qe ([0,l- J R]U(uf =1 [2j-(l-i?),2j + (l-i?)])U[2iV-(l-i?),2Ar])x[-l,l]. 

Note that if we evolve the system starting from an initial state containing a trapped 
particle along any sample path, the system is going to contain a trapped particle at all 
times. 
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Let St C fl be the set of all states with trapped particles. 

Theorem 1. Suppose there exists a probability measure \i invariant under <P T with 
I-i(St) = 0. Then fi is absolutely continuous with respect to the Lebesgue measure 
m and ergodic. 

Remarkl. Since obviously n(f2k) 7^ for all k, our assertion that fi is absolutely 
continuous with respect to m means that on each Q^, it has a density with respect to 
m k . 

Corollary 1. If v is an ergodic measure for <P T , then for some k > 0, v is supported 
on Uj°_ k f2j and can be represented as a direct product of two measures, v = /i x tt, 
such that fx is the unique absolutely continuous ergodic measure with /J.(St) = as 
in Theorem\l\and tt is a singular measure supported on the states in ilk traced by the 
trajectories ofk trapped particles. 

Equilibrium Case: 

Suppose the system is coupled to two equal heat baths, characterized by temperature 
T and injection rate g, i.e. the injection process at each bath is Poisson with rate g, the 
distributions for the positions of injections are uniform on an 7^, |7l| = |7^| = 
I7I, and upon injection, a particle is assigned a random velocity v sampled from the 
distribution 

ce- m ^\'\v\coe(<p)dv, c= gi^g ; 
where f3 = l/T, m is a particle's mass, and ip E (— f , 

Theorem 2. The invariant probability measure fi characterized by the properties below 
is ergodic. 

(i) The number of particles in the cell is a Poisson random variable with mean 

-area(r) gy/m 
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(ii) fi has conditional densities Ck<Tkdmi~ on flk where Ck is a normalizing constant 
and for X G flk, 

a k (X) = e-^ £f =1 4+mi:U ^#). 

Proof. For 77 = = 1, the invariance of /i is shown in Q. The argument can be 
generalized for any < 77 < 00. Since fi <C m and /i(S't) = 0, ergodicity of // follows 
from Theorem Q] □ 

The remaining part of the paper is devoted to the proof of Theorem Q] In order to 
simplify the exposition, we chose to present the proofs of all the technical lemmas for 
the situation 77 — 1 only; all proofs could be generalized for any < ?/ < 00. Aside 
from technical lemmas, the argument deals with any 77, < 77 < 00. 
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3. Outline of Proof of Theorem [T] 

In this section we state three propositions and show that they imply Theorem [1] Our 
proof uses some ideas from jT|. 

Idea of proof: Proposition Q] states that absolutely continuous measures stay abso- 
lutely continuous under <P T . Propositions[2]and[3]imply that any singular measure even- 
tually acquires an absolutely continuous component when evolved under <P T . If follows 
that the singular component of any invariant probability measure with [i(St) = must 
be zero. Propositions |2]and|3]also imply that m-almost all initial states must belong to 
the same ergodic component; ergodicity follows by the absolute continuity of fx. 

Denote by 

- Ssc the se t of all states in fl such that a simultaneous collision with same disk 
occurs under the evolution of the system with no particle injections, i.e. for some 
t > 0, q\ € dDk and € dDk for i ^ j. Note that the evolution of the system is 
not defined after such a collision. 

- Stc the set of all states in SI such that a particle stops under the evolution of the 
system with no particle injections, i.e. 3t such that Vr > t, v* = 0. When r\ = 1, 
this situation occurs when a particle hits a stopped disk tangentially; when 77 ^ 1, it 
occurs when a particle hits a disk with angular velocity ui = ^ T |~|_\^ t tangentially, 
where is a tangential component of the particle's velocity upon collision. 

Let S = S T U S S c U S T c- 

Definition 5. A state lefi is called admissible ifX^S. 

Given t > 0, the probability that no particles are injected on time interval [0, t] is 
positive. If \i is an invariant probability measure with /i(Sr) = 0, then it cannot give 
positive measure to either Ssc or Stc, i- e - n(S) = (i(St U Ssc U Stc) = 0. 

It follows that if we start with any measure v <C /1, the push forward of v under the 
Markov process <P T , (# T )*z/, is well defined for all r > 0. 

For any measure v, denote by and v±_ the absolutely continuous and singular 
components of v with respect to the Lebesgue measure to, i.e. v = v<^ + v±, <§C to 
and u± _L to. 

Proposition 1. If v ^ m, then (f? t )*y is well defined for any t > and (<P t )*i / *C m. 
In particular, (^t)*(/ i <) ^ m for any t > 0. 

Consider a sequence of particle injections c = (ci, • • ■ , c„) at times < ti < £2 < 
• • • < t n < T such that at time tu a particle enters the system at location £j 6 7^, 
at angle 5j € (— f, f), and with speed s» € (0, 00). Then, assuming no simultaneous 
collisions with the same disks occur, one can generate a sample path a defined on [0, T] 
in which one starts from state X E Q and injects particles into the system according to 
c. 

We call C a canonical neighborhood of c if there are disjoint open neighborhoods 
Tj of ti contained in [0, T], Si of £j, A; of 6i, and Sj of such that for each sequence 
of injections d E C, exactly one particle is injected in each Tj, with position in Qi, 
angle in Ai, and speed in Si. No other injections occur in the time interval [0, T], 

We call £ a canonical neighborhood of a if there exist an open neighborhood U 
of X and a canonical neighborhood C of c such that each sample path in starts 
with an initial condition in U, is generated by a sequence of injections from C, and no 



■s 



Tatiana Yarmola 



simultaneous collisions with same disks occur. Note that for any Y £ U, if y is a set of 
all sample paths on [0, T] starting at Y, then the probability that a sample path from y 
belongs to £ is positive. 

Propositions 12 and |3] split the problem of acquiring density for singular measures in 
the following way: Proposition|2]deals with the simplest situation of acquiring density 
for a point measure supported on an particle-less initial state Yq £ S7q\ Proposition [3] 
provides a sample path from any admissible initial state to a particle-less state, from 
which a density can be acquired using Proposition [2] 

Proposition 2. Given a state Yq £ Qq with a condition that if rj — 1, all disks have 
nonzero angular velocities, there exist an open neighborhood Uq o/Yq, time Tq, and a 
set Aq C J?o with itiqIAq) > 0, such that for any Y £ Uq, [(^t )*<^y']< has strictly 
positive density on Aq. In particular, [(^r )* < ^y]<c( 1 ^) 7^ 0- 

Proposition 3. Given an admissible state X £ Q, a state Yq £ Qq, and a neighborhood 
Uq C fio ofYo, there exist time T, a sample path a on [0, T] that starts at X and ends 
at Yq, and a canonical neighborhood S of a, such that each sample path in S ends in 
Uq. 

Proof ofTheorem\l\assuming Propositions\J\\3\ and\2\ 

Propositions |2] and |3]imply that for any admissible state X £ Q, 3 a neighborhood 
U of X, such that V Y £ U, [(<Pt+t )*<V]< has strictly positive density on Aq and, in 
particular, [($t+t )*5y]<(O) ^ 0. 

Assume fi±(f2) ^ 0. Since \i is invariant with \i{St) = 0, fi±(S) = H>(S) = 0. 
Therefore [(<I>T+To)*f 1 ±-]<g:(ft) 7^ 0. Applying Proposition Q] we conclude that Vt > 
T + T ,[(#t)./i±]<(rt)^Q. 

Clearly (# t )»/i = [(#*)*(//<)]< + [(**)• (m«)]x + [(*«)* + [(**)* (mi)]-L- 

By Proposition!!] Vt > 0, [(^t)*(/i<)]±(^) = 0. Therefore, for t > T + T , 

[(* t )./*]<(l2) > M<(^), 

which contradicts the invariance of \i. This proves the absolute continuity of /j, with 
respect to the Lebesgue measure m. 

Assume /_ti and fi2 are ergodic measures with hi(St) — /^(Sr) = 0. Then /ij <C 
to, i = 1,2. Suppose there exists a Borel function cp such that ci = f n (pdfii ^ 
J n (fid/12 = C2- Then by the Random Ergodic Theorem, there exist A4 G Q, i = 1, 2, 
with jUi(vlj) = 1 such that for every X £ A4 the ergodic averages are equal to Cj for a.e. 
sample path starting from X; and m(Aj) > 0, since fa -C to. Since for any admissible 
state X £ f2, [(<Pt+t )*Sx}<^ has strictly positive density on Aq, the random ergodic 
averages for all the admissible states are equal for measure 1 sets of sample paths. Since 
771(6*) = 0, either m(Ai) — or 771(^2) = 0, a contradiction. Therefore fit = H2 = /■*• 
□ 

We prove Proposition [TJin section|4] Proposition |2]in section|5] and Proposition[3]in 
sections [6] and [7] 

4. Proof of Proposition [T] 

In order to show that, given v <C m and t > 0, ($t)* v is defined and ($t)* v "C 771, 
we would like to consider a countable number of subcases depending on how many 
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particles entered on time interval (0, t] and show that for each subcase the time-t push 
forward of v is well defined and absolutely continuous. 

Suppose n particles enter on time interval (0, t\. For n > 1, let C n = [(0, t] x (7^ U 
7fl) x ( — § 1 § ) x (0, oo)] n / ~ be the set of all possible injection parameters for these n 
particles modulo permutations and let p n be the natural probability measure on C„, i.e. 
the product of appropriate injection distributions modulo permutations (see subsection 
11.2b . Given t\, ■ ■ ■ , t n -\ with = to < ti < • • • < t n = t, denote by Ct lt ... ,t„_! the 
subset of C n such that j th particle is injected on (tj-i,tj]. 

For any r > 0, denote by (^™) >t i/ the time-r push forward of v under the random 
dynamics assuming that n particles entered on time interval (0, r], if defined. Given a 
subset C C C n of injection parameters, denote by (^" ,c )*j/ the time-r push forward of 
v under the random dynamics assuming that n particles entered on time interval (0, r] 
with injection parameters from C, if defined. 

Lemma 1. If v -C to, (<P^)^iy is defined and absolutely continuous with respect to m. 
Moreover, for any n > 1 and a choice of 't±, • • • , f n _i with = to < t\ < ■ ■ ■ < t n — t, 

(<I> t ' 41 ' )* v is defined and absolutely continuous with respect to m. 

Proof of Proposition\J\assuming Lemma\J\ 

Let B = {X € i? : P l x is not defined }. Then for each X 6 B, there exists a positive 
measure set Cx of injections such that, if we start at X and follow any sequence of 
injections c G Cx, the dynamics is not defined up to time t, i.e. a simultaneous collision 
with same disk occurs before time t for each such injection. 

Assume m(B) > 0. Define A = {(X, c) : X e B; c <E C x }- Then by Fubini's the- 
orem (to x p n )(A) > 0. Let (X, c) be a Lebesgue density point of A with X admissible 
and c having no simultaneous injections ((X, c) exists since 771(5) = and probability 
of simultaneous injections is zero). If t ei , ■ ■ ■ , t Cn are the injection times for c, choose 
any t\, ■ ■ ■ , t n -i such that = to < t Cl < t\ < ■ ■ ■ < t„_i < t c „ < t n =t (c cannot 
be an empty sequence of injections since is defined by LemmaQ]). Let U be any 

neighborhood of X such that any Y € U is admissible (the set of admissible states is 
open and dense in ft). Then (to x p n )[A f~l {U X C^,... t n _ x )] > 0, while by LemmaQ] 
it must be zero. A contradiction. Therefore m(B) — and is defined. 

Suppose (#™ )* 1/ has a nonzero singular component supported on a set F C Q. Then 
m(F) — and to{X g i7 : P x (F) > 0} > 0. By an argument similar to the above we 
get a contradiction to LemmaQ] Therefore (#™)*f is defined and absolutely continuous 
with respect to to, which completes the proof of PropositionQ] □ 

Proof ofLemma\l\ 

Since Ct lt ... I t n _ 1 C C n is such that one particle is injected on each (tj-i, tj) and 
the probability of two particles entering at the same time is zero, to prove LemmaQ]it is 
enough to treat two situations: particles are injected on (0, t] and 1 particle is injected 
on(0,t]. 

Case 1(0 particles injected). Since m(S) = 0, where S is a set of non-admissible 
states, and v <C to, the time-r push forward of v is well defined in this situation. 

The system behaves as a closed one until some particles exit. We would like to split 
the possible situations according to the number of particles that exit: 

0: no particles exit. 

Let Aq C fl be the set of states such that no particles exit on time interval (0, t]. 
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Decompose Aq = UAq such that A§ = Aq n Qk- Then v\ A k <C m/- and must stay 
so at all times (0, t], i.e. (<&®)*(v\ A k ) -C Vr £ (0, t]. Indeed, if for some t < t 
andaBorel set B, [($°)*(i;\ A k)] ± (B) ^Obutm fe (B) = 0, then (v\ Ak )($_ T B) > 
[(^°)*(i>\ A k)]±(B) ^ 0, while m k {^- T B) = m k (B) = 0, a contradiction to the 
absolute continuity of v\ A k. Thus {<P^)^{v\ Aa ) <C m. 

1: 1 particle exits or several particles exit at the same time. 

Let A\ be the set of states such that 1 particle exits on the time interval (r, t] and 
does not collide with dr \ (7^ U 77?) on that time interval. If several particles 
exit at the same time, let A\ be the set of states such that neither of these parti- 
cles collides with dr \ (7^ U 7^) on the time interval (r, t]. Then until time r, the 
dynamics is equivalent to the dynamics of the closed system and {<P®)*(v\ A t) <C 
m. Since after the time r the particle(s) do not collide with the disk, the parti- 
cle^) coordinates are independent from the rest of the system and (<P^ )*(^U^) = 
(^j_ t )*(P[(<?°)*(^| j 4t)]) <c m, where P denotes the projection of the measure 
on the remaining coordinates in the system. Since the statement is true for any 
re [0,t), (*?)*(i/Ui)< m, where A 1 = U TeM Al. 

n: particles exit at n different times with possibly several exiting at each time. 

Let A n be the set of states such that n particle exits occur on the time interval (0, t], 
with possibly several particles exiting at the same time. Then 

A n — Uo={ <t 1 <...<t ra =tAt li ... 

where A tl .... ,t n _ 1 denotes the set of states such that j th particle exists on the time 
interval (tj-i,tj], 1 < j < n. Each A tl: ... i t n _ 1 can be treated as in the previous 
situation by applying the argument consequently to each time interval (tj-i,tj], 
1 < j < n. Thus (#()*Ha„) < m - 

We conclude that (<l^)^v <C m. 

Case 2(1 particle injected). In order treat this case, we need to study the injection 
process in detail. We would like to start with a simplified situation when the injected 
particle does not collide with disks or walls on time interval (0, t]. We start by studying 
the injection process for the cell r occupying the x > half plane of xy-plane bounded 
by the vertical wall along the y-axis with an opening (—7/2, 7/2) along the y-axis. 

Lemma 2 (Particles are injected with 4-dimensional uncertainty). Let P be as above. 
Assume a particle is injected through the opening at a random entrance time r G 
(0, 00), with a random position £ G (7/2, 7/2), at a random angle 5 G (— 7r/2, it/2) 
and with a random speed s G (0, 00). Assume also that the distributions for entrance 
time, position, angle and speed are finite and positive on all the given intervals. Then 
the measure that describes the probability of finding the particle at certain position 
(a;*, ?/*) with certain velocity (w* , u* ) at any time t > is absolutely continuous with 

respect to the Lebesgue measure on r xl 2 and has positive density everywhere. 

Proof of Lemma\2\ The lemma holds if the mapping /' : (0, 00) x (7/2,7/2) x 
(— 7r/2, 7r/2)x (0, 00) -> (0, 00) x (— 00, 00) x (0, 00) x (— 00, 00) such that (r, £, S, s) -> 
(x*, ?/*, i>* , Vy ) is a diffeomorphism. This boils down to verifying that the Jacobian de- 
terminant is nonzero everywhere and that the map is surjective. Simple computations 
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yield that: x* = (t — t)scos(S), y* = — (t — t)s sin((5) + £, = scos(<5), and 
Vy = —s sin(<5) and thus the Jacobian matrix is: 

/ -s cos(<5) -(t-r)s sin(5) (i-r)cos(5) \ 
ssin(<5) 1 — (t — t)s cos((5) — (t — r) sin(j) I 
-ssm(6) cos(<5) 
V -scos((5) -sm(S) J 

and its determinant is — s 2 cos(#), which is nonzero on the domain of definition of /*. 
Given (a;*, y*, v\ , let 

s = ^(4) 2 + K) 2 . 5 = arctan(-»*/i£), 

x* sin(<5) + y* cos((5) y* — £ + ts sin(<5) 

cos(<5) s sm(6) 

Such (r, £, (5, s) clearly maps to (x*, y*, u* , «*). □ 

In Lemma |2] the domain is such that an injected particle cannot collide with the 
boundary of r at any time. In the real situation, when a particle is injected into the 
playground r, it might collide with dr in an arbitrarily short time, depending on the 
injection parameters. 

Let Cf be the set of all particle injection parameters (r, £ , 5, s) such that one particle 
enters on time interval (0, «], < K < t, and does not collide with dr on time interval 
(0, re]; no other particles are injected on time interval (0, t]. Since we assumed that all 
injection distributions have positive density everywhere and Cf clearly has a nonempty 
interior, by Lemma [2] the measure that describes the particle location and velocity at 
time re is absolutely continuous with respect to the Lebesgue measure on T x I 2 ; 
denote this measure by ire* ■ 

Suppose we start with a measure v <C m and assume that one particle is injected 
on (0, t] with injection parameters drawn from C\. Then until time re the coordinates 
associated with the injected particle are uncoupled from the coordinates of the rest of 
the system. Thus v evolves as in Case 1, is defined, and ((p^^v <C m. 

Therefore 

(#i' cr ),f - [ircn] x [(<£°M « m. 

The injected particle have entered by time re and no other particles enter on time inter- 
val (0, t]. Thus on time interval (re, t] the system behaves as in Case 1, implying that 

(^J' 1 )*v <C to. Since the above result is true for any re € (0, t], {p\) t v <C to. □ 
5. Proof of Proposition [2] 

We are going to "acquire density" starting with an initial state Yq £ J?o by hitting each 
disk once with a particle possessing a "4-dimensional uncertainty". Upon each such 
collision, each disk would acquire a "2-dimensional uncertainty." That is, we will show 
that there exists an open set of injections C of N particles, such that the j th particle 
hits disk Dj once and exits the system with no additional collisions by some uniformly 

selected time To. Moreover, (^^' C )*^y Q <C to. It will turn out that same statement is 

true for nearby states Y £ Hq with (<P^, g C ) :t SY varying continuously with Y. 
In order to prove Proposition |2] we need the following lemmas: 
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Lemma 3 (Continuity). Let B C r x R 2 be an open ball of positions and velocities 
such that for any (x,y,v x ,v y ) £ B a particle with position (x, y) and velocity (v x , v y ) 
at time t = it is going to hit disk Dj within time tq while not meeting dT again on 
the time interval [0, t{\, t\ > tq. Note that the collision with Dj must be non-tangential 
since B is open. Let Fj be an open ball of angular positions and velocities of disk Dj 
at time 0. For any t £ (r , n), define f t :BxFj^rxR 2 x dDj x M such that 
ft{x, y, v x , v y , <p, uj) = (x , y , v x , v y , ip ,uj ) gives the positions and the velocities of 
the particle and disk Dj at time t. Then f t is continuous for any t £ (to, ti). 

Lemma 4 (Acquiring density for a disk). Let B C r x R 2 be as in Lemma\3\ Assume 
at time disk Dj has position <p and angular velocity to. Define Hff' M : B — > S 1 x R, 
(x, y, v x , V y ) — > ((p , L) ) to be the mapping of the particle position and velocity at 
time zero to the disk position and velocity at time t £ (tq, t\). Let v be a measure on 
B equivalent to the Lebesgue. Then the push forward measure {^f^)^^ is absolutely 
continuous with respect to the Lebesgue measure on S 1 x R and has positive density 
on some open set. 

Proof of Proposition\2\assuming Lemmas\3\and\4\ 

Denote the angular positions and angular velocities of the disks in state Yq by 
(ipi,u!i), ■ ■ ■ , (lpn,ujn). By assumption, if T] = 1, uij ^ 0, 1 < j < N. 

Suppose we inject a particle at time r > with some initial position £ £ jl, initial 
angle tp £ (— \ , and speed s £ (0, oo) arranged in such a way that it first collides 
non-tangentially with Dj at the top, (2j — 1, R), and exits the system with no additional 
collisions; this is possible because ojj ^ 0, 1 < j < N, if r\ = 1. Since the collision 
with disk in non-tangential, by Lemmas [2] and [3] there exist open neighborhoods Ij of 
ujj and Cj of (r, £, ip, s) such that for each ui' £ Ij and (t, £, ip, s)' £ Cj, the injected 
particle follows a nearby path in r, hits Dj once, and exits the system with no additional 
collisions within some time Tj . 

Suppose we subsequently hit each disk with a particle as described above. Define a 
neighborhood U of Y by U — dD\ x • ■ ■ x dDjy x h x • • ■ x I N . 

Let C = nf =l Cj and T = J2j=i Tj- Then by Lemmas and gj for any Y £ U , 

(<& T ' )*5y *C tuq. Here (^ T ' ) x Sy denotes the time-T push forward of 5y, provided 
that exactly N particles enter on time interval (0, T] allowing only injections with pa- 
rameters in C. Let Ay be the set of states on which (<? T : )*8y has strictly positive 
density. By Lemma|U each Ay contains an open set and by Lemma[3] Ay vary contin- 
uously with Y. Therefore there exists an open neighborhood Uq of Yq, Uq C U, such 
that Aq = r)y<zu Ay contains an open set; clearly mo(Ao) > 0. This completes the 
proof of Proposition |2] □ 

Proof of Lemmas\3\and^\(t] = 1). 

Assume the disk is of radius R and is centered at (0, 0) coordinate of the xy-plane 
and the particle and the disk have initial coordinates (x, y, v x ,v y ) and (tp, uj) respec- 
tively. 

Denote by r the collision time with the disk, by 9 the angular position of the collision 
point on the disk measured counterclockwise from the positive direction of the x-axis, 
and by v t and v± the tangential and the normal velocities of the particle upon collision, 
with v± representing the velocity after the collision, i.e. pointing outwards. Then, 

Vt = —v y cos(0) + v x sin(0) 

vj_ = —v y sin(#) — v x cos(#) 
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Rcos(9) — x Rsin(6) — y 



and 

One can rewrite 



2 i 2 
V x + V y 



Rv t = v x y - v y x 



Rvx = -t(v x + v ) - v y y - v x x = J v% + v. 



Following the collision, at time t E (to, ti), the coordinates for position and velocity 
of the particle are: 

u* = Vj_(tv x + x)/R + lo(tv v + y) 
v l = vj_(rv y +y)/R- lo{tv x +x) 

x l = tv x + x + v' x (t - t) 

y l = Tv y + y + Vy(t - t) 

The angular velocity of the disk at time t is w l = Vt/R and the angular position is 
= [ip + lot + Vt(t — r) / R] mod 2irR. Clearly / t is continuous, which completes 

the proof of Lemma[3] 

In order to prove Lemma |4] we would like to study the matrix of the derivatives of 

S 1 xR, (i,i/ l ^^)->(^w t ),t6 (ro,ri).Itis: 



/ dip 1 diJ_\ 

dx dx ' 

3y* ckJ_ 

dy dy 

dy- ckJ_ 

, V a/ , 

\ dv y dv y ) 



+ w\* T > 



^y_( i ._ rr \ _^v_\ 



X 



V 2 -\-V 2 

(-?f-2%T-y)(w-^) 



R? (* 



V 



vl+vl 



R 2 (* T ) ft 2 / 



where T denotes the transpose. 

Suppose the two rows of the derivative matrix (i.e. columns of the non-transposed 
matrix) are linearly dependent, then there exists a constant c such that: 



R 2 (- 



VtVy 
V± 



v x )(u 







-\- IP" 



(t 



v 4- v 



+ (*-r) 



i? 2(^_2^r- a ;)( W -f) 



-I- 



Thus either (_Rcj — ft) = or 



-(*- T ) = " + ( t - T ) 



V 2 4- V 
u x ' u y 



Vt | Vx 

V±_ Vy 



Vt 
V± 



vt _ 2v x t _ x _ vt_ 2v y T y 
v± y y v± xx 
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From our assumptions it follows that x, y, R, v x ,v y ^ 0; so the derivative matrix has 
rank 2 unless Rcu = vt- 

Let A u = {(x,y,v x ,v y ) : Ruj = v t = j}(v x y - v y x)}. Clearly v{AJ) = 0. 
For each point (x, y, v x , v y ) G B \ A^, the derivative matrix has rank 2. Therefore 
the push forward of v under W^'" must be absolutely continuous with respect to the 
Lebesgue measure on S 1 x K and has positive density on an open set \Pf' u {B \ A^). 
This completes the proof of Lemma|4] □ 

6. Flushing Particles Out 

This section is the I s * step in the proof of Proposition|3] Here we will show the follow- 
ing: 

Proposition 4. For any admissible state X, there exists a sample path ax that starts at 
X and ends at some Xq € J?o 

In order to drive the system from state X with possibly many particles to some 
particle-less state Xq, we have to ensure that each particle in X traces a path in r from 
its initial position to one of the exits. Following the ideas in |4), we will describe a class 
of projected particle paths traced in _T and show that each can be followed provided 
that disks have appropriate angular velocities upon collisions. Then we will establish 
that by injecting particles with appropriate initial conditions we can change the angular 
velocity of any disk to any given value in an arbitrarily short time. That will enable us 
to force a particle along a projected particle path by setting the angular velocities of the 
disks to appropriate values before collisions. 

Definition 6. A proper projected particle path is a continuous curve 7 : [0, 1] — > r, 

s 1 y j(s), such that 

1. 7 consists of a finite sequence of straight segments meeting at dF. 

2. The incoming and outgoing angles of two consecutive segments of 7 meeting OTq 
are equal. 

3. Only 7(0) and 7(1) can be in the openings and 7^. 

4. 7 /s nowhere tangent to boundaries of the disks UjL 1 dDj. 

Remark 2. Note that a proper projected particle path is allowed to have any non-tangential 
'reflections' off the boundaries of the disks \J^ =1 dDj. An example of a proper projected 
particle path is shown in FigQ] 

Lemma 5 (Existence of a proper projected particle path). There exits a proper pro- 
jected particle path from any point (x, y) £ dDj to one of the exits Jl or jr. 

The statement in Lemma [5] is rather obvious for the geometry we consider. We in- 
clude a proof for completeness purposes. 

Proof of Lemma\5\ Let S be such that < 5 < 2R ^ F , R . If y-coordinate of (x, y) is 
not in [—5, S], i.e. \y\ > 6, then by a simple geometric argument there exists a proper 
projected particle path from (x, y) to the top of the disk Dj, (2j — 1,R) (or the bottom 
(2j — 1, —R)), which makes several collisions with the upper wall [0, N] x {1} (or the 
lower wall [0, N] x {—1}). By appending to the proper projected particle path above a 
segment in -T that connects (2j — 1, ±fc) to either jl or 7fl. > we prove Lemma [5] for 
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the situation when \y\ > S. If y E [—5, 6], then there exits a segment in r that connects 
(x, y) to the point on the appropriate nearby disk or exit with y-coordinate between 5 

and 2R 2-r R ■ Indeed, by the system's geometry, a tangent line from (2j— 1±R, 0) to the 
appropriate nearby disk (if applicable) intersects the disk with ^-coordinate ± 2R ^ / _} F! — ■ 

a 

In order to force a particle to follow a proper projected particle path, we have to 
ensure that upon each (non-tangential) collision with a disk, the disk has appropriate 
angular velocity: the total velocity of the particle after collision is the vector sum of 
v'j_ = —vj_ and v' t = Vt — Y^{ v t — R<^>) and must be parallel to the next segment of 
the proper projected particle path. We would like to establish that by injecting particles 
with appropriate initial conditions we can change the angular velocity of any disk to 
any given value in an arbitrarily short time. 

Lemma 6 (Controlling angular velocities of disks in arbitrarily short times). Sup- 
pose disk Dj rotates with angular velocity w and none of the particles inside the system 
will collide with any disk before time r > 0. Given any ui' r there exists a sequence of 
particle injections on time interval (0, r) from the left bath such that: 

— at time t the disk Dj has angular velocity uf ' , 

— at time r all the injected particles have left the system, and 

— on time interval (0, r) the injected particles follow admissible paths and only hit 
disks D\, • • • , -Dj-i with the exception of one collision of one particle with disk Dj. 

The same holds for the right bath with appropriate disk renumbering. 
We prove the Lemma|6]in subsection !6.3l 



6.1. Proof of Proposition^ no tangential collisions. By definition of an admissible 
state, one of the following holds for each particle in X under the evolution of the system 
with no particle injections: 

1. there exists a finite time t > such that the particle exits the system at time t 
and does not collide with any disks on time interval [0,t], i.e. q f e Jl U 7_r and 
<f ujLx^i Vt e [0,t]; 

2. there exists a finite time t > such that at time t the particle collides with a disk 
non-tangentially and no other disk collisions occur on time interval [0, t], i.e. <f g 
uf =1 dDj with v]_ ^ and q T <£ Uf =1 dDj Vr e [0, t]\ 

3. there exists a finite time t > such that at time t the particle collides with a disk 
tangentially with v x ^ and no other disk collisions occur on time interval [0,t], 
i.e. q l G \Jf =1 dDj with v]_ = and v x ^ and q T uf =l dDj Vt G [0, t]. 

In this subsection we are going to assume that for each particle in X either 1 or 2 
holds, i.e. under the evolution of the system with no particle injections each particle 
in X either exits the system or collides with a disk non-tangentially. We will treat the 
situation with tangential collisions in subsection !6.2l 

Suppose X contains k particles and the above assumption is satisfied. Then using 
Lemma|5]for the j th particle in X, 1 < j < k, we can assign a proper projected particle 
path 7j from the particle's initial position Vj) to one of the exits. If we force the 
particle to follow this path, the times of all collisions are fixed. Indeed, the unique angu- 
lar velocity of the disk keeps the particle on the path for each (non-tangential) collision, 
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implying that the speeds with which the particle traces segments of the path and there- 
fore the collision times are uniquely determined from the path. If j th particle follows 
jj, let rf , • • • , T^tj) ^ e times of collisions with disks, D^(j t x) > ' " " > ^k{j,n(j)) an< ^ 
let ojj , • • ■ , <^/j) be the required angular velocities. 

Assume first that all rf (for all particles in X and all collisions) are different. Then 
direct application of the Lemma [6] between collisions guarantees the existence of a 
sample path ax from state X to some state Xq e S7q. 

If some t\ happen to coincide, simultaneous collisions with same disks might occur, 
making us unable to construct ax with our choice of paths. We would like to show that 
we can always choose a collection of nearby paths from the particles' initial positions 
to the exits such that no simultaneous collisions with same disks occur. Note that we 
cannot, in general, avoid all simultaneous disk collisions since we have no control over 
the times of the first disk collisions of the particles in X. First collisions, however, 
cannot occur simultaneously with same disks since X is admissible. 

We treat the possibility of simultaneous collisions with different disks as follows: If 
several, say n, collisions are about to occur in time r after the previous collision time, 
all with different disks, we can set the angular velocities of the disks by the Lemma [6] 
in the order of decreasing disk index, in a fraction of time r, ^, each: the Lemma [6] 
guarantees that the disks with indexes larger than the one we set the angular velocity 
for are left untouched. 

In the remaining part of the proof we show that we can choose a collection of paths 
for particles in X from their initial positions to the exits such that no simultaneous 
collisions occur. We start with a description of a set of paths to choose from for each 
particle in X. In the following Lemma, we assume that at time the j th particle is the 
only particle in the system in order to ensure that the system is defined at all times. 

Lemma 7. Assume that under the evolution of the system with no particle injections 
the j particle collides with a disk non-tangentially in some finite time t > 0. Then 
there exist open neighborhoods If ofu>{, ■ ■ ■ , I n ^j\ °f^n(j) suc ^ that for any choice of 
angular velocities (oj^)' € I(, ■ ■ ■ , (^yp' € P n ijy if a particle starts with its initial 
position (q J ,v J ), it collides with D 3 k{jl y ■ ■ ■ iD 3 k{jn{j)) set to (lu{)' , ■ ■ ■, {u 3 n{j) )' , and 
exits the system. 

Lemma|7]follows from Lemma[3]and the continuity of the billiard flow. 

Lemma [7] guarantees that the j th particle exits the system for any choice of angular 
velocities from • • • , P,^ If a particle starts at (q^v^), for different choices of an- 
gular velocities in I{ for disk Df-y^, the particle collides with disk D^j^) at different 
times. In fact, one gets an open set of possible collision times since the set of possi- 
ble positions of the particle at any time r from the collision with disk D k ^ ^ forms a 
broken line: before any collisions, the particle's positions are 

tv 1 = t(v' x + v' t ) = r(v' x + [v t - - R(<4)% = 

= t(v' ± + [ Vt - ^-{v t - M)]) + - (<4n 

1+7] 1+7] 

where (uj{)' varies through I\ ; and upon reflections from straight walls the straight line 
of positions becomes a broken line. If we fix a specific (ui{y € I{, similar argument 
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shows that different (lu^Y G -T| f° r disk Dk(j,2) yield a range of collision times with 
D 3 k Q 3 \! and so on. 

Therefore, can always pick a collection of (u> 3 )' G such tnat no simultaneous 
collisions occur for all particles in the system. Given such a choice, each particle in 
X follows a path from its initial position to an exit. Setting the angular velocities of 
the disks to (uj 3 )' at appropriate times by following appropriate sequence of injections 
provided by the Lemma [6] we construct a sample path ax from state X to some state 
Xq in [2q. This completes the proof of Proposition [4] for the initial states X in which 
particles either hit a disk non-tangentially or exit the system. □ 



6.2. Tangential Collisions. In this section we will prove Proposition |4]for the general 
case, i.e. when for some particles in X the first collision with a disk might be tangential. 
Since X is admissible, such collisions must occur with v x ^ and without simultane- 
ous collisions with same disks. 

Suppose the first collision of the j th particle is tangential with v x ^ 0. By the 
Lemma|6]we can ensure that, upon collision, R times the angular velocity of the disk is 
equal to the velocity of the j th particle, Ru = v = U{. After such a collision the particle 
continues along the straight line with the same velocity as if no collision has occurred. 
From that point, the particle will either exit the system, hit a disk non-tangentially, or 
hit a disk tangentially again. In the third situation we would like to set the angular 
velocity of the disk to be equal to the velocity of the particle again and continue the 
process. Since v x ^ and |t>| are kept constant under subsequent iterations of the 
third situation and upon collisions with straight walls, the third situation can occur at 
most finite number of times; then either the particle will either exit the system or will 
collide with a disk non-tangentially. In the later case, by Lemma[5] there exists a proper 
projected particle path from the non-tangential collision point to an exit. Denote the 
path traced r by the j th particle in this construction by jj. 

Again, the times of all collisions of j th particle along jj are fixed; let , ■ ■ ■ , t^, -s 
be the times of collisions with disks, D k (j^ , ■ • • , Dk(j,n(j)) an d let • • • , w^y) be 
the required angular velocities (with Ruif = v = v t at tangential collisions). In the 
following Lemma, we again assume that at time 0, the j th particle is the only particle 
in the system in order to ensure that the system is defined at all times. 

Lemma 8. Assume the j particle has m(j) > tangential collisions before it has 
a non-tangential collision or exits the system. Then there exist open neighborhoods P 1 
of ■ ■ ■ , I n (j} °f ^lid) suc h that for any choice of angular velocities (cj{) G I{, 

■ ■ ■ , (uj° n ^Y G I^jy if a particle starts with its initial position and velocity (qj,Vj), it 
possibly collides with ffcy,!) , • • • , -Dfc(j,m(j)) set to (w^)', ■ ■ • , {uj^^Y , collides with 
D-j i . , D, , ■ ,.„ set to (aj J ,-,,,)', • • • , {uj 3 , s )', and exits the system. 

Proof. By Lemma|7]there exist open neighborhoods I 3 m(j)+1Q of uj 3 m{j)+1 , • • • , I 3 n{j) 
of uj 3 , .s such that for any choice of angular velocities (u> 3 ,••,,,)' G II , ■■, , , , • • ■ , 
{u} 3 n ^Y G ^n(j) if a P^tide starts with (qj,Vj), it possibly collides with D k ^ ^, 

■ •• . D k{ 3 , m { 3 )) set to uj{,--- , w^yj, collides with D k{j ^ m(j)+1) , D k{jMj)) set 
to (to 3 , s , , )' G II i -s , , n , • • • , (to 3 , .> )' G P i .v n , and exits the system. 
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Fig. 2. uj' = case: the incoming trajectory L is marked by a solid line; outgoing V - by a dashed line; 
velocities are labeled after collisions 

At tangential collisions, we chose to set R times the angular velocities of the disks to 
be equal to the velocities of the colliding particles, Rlj = v = vt-lf we vary the angular 
velocities of the disks at tangential collisions keeping sign the same, the particle still 
follows the same path and exit the system. □ 

If a particle starts at (q 3 0l Vq) and the collision with disk Dk(ji) is tangential, letting 
the angular velocity of £ ) fe(j l i) to vary through If, we get an open set of possible col- 
lision times with D^j^) since the particle follows the same path, only with different 
speeds. The remaining tangential collisions are treated analogously and non-tangential 
collisions as in subsection l6.ll 

Therefore we can always pick a collection of (ujf)' G I\ such that no simultaneous 
collisions with same disks occur for all particles in the system. Given such a choice, 
each particle in X follows a path from its initial position to an exit. Setting the angular 
velocities of the disks to (cuf)' at appropriate times using the Lemma [6] we construct 
a sample path ax from state X to some state Xq in J7o.This completes the proof of 
Proposition^ □ 

6.3. Proof of the Lemma\6\(rj = 1). The proof is by induction on j, 1 < j < N. 

j=l: There are many ways to treat this case, but we choose a very specific one that 
will be a useful step in treating the induction step. Our method here is identical to the 
one in proof of Lemma 4.3 in J4). 

Assume first w' — 0, i.e. we want to hit D\ radially. Send a particle parallel to 
the x-axis to hit the disk D\ at (1, 0). If the initial velocity v = (v x , 0) is big enough 
compared to Ruj and the distance from jl to D\, the particle will be able to exit the 
playground without hitting dT \ (7^ U 7^) again. Note that the larger v x is, the smaller 
is the angle of reflection; we can introduce any bound on the angle by choosing v x large 
enough. When w' 7^ 0, we send a particle such that it also hits the disk at (1, 0). We 
can introduce any bound on the angle of incidence in a similar way. Clearly, by making 
v x sufficiently large and large enough compared to Rw', we can complete the above 
procedure in an arbitrarily short time. 

Induction step: Assume the lemma holds for all j < k. We would like to show that 
it also holds for j = k + 1. We want to send a particle with velocity v = (v x , v y ) such 
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that it first hits disk Dk at (2k — 1, —R) without hitting dT along the way. The velocity 
after collision is (RuJk, —v y ), where ui k is the angular velocity of disk Dk, which we 
can set to any value in arbitrarily short time by the induction assumption. 

Consider first the case to' = illustrated in Figure|2] Then we want to hit disk Dk+i 
radially by following the unique trajectory that reflects from the lower boundary of the 
playground r, [0, 27V] x {—1}, once before hitting D^+i; call this trajectory L. If is 
the angular position of the collision counting counterclockwise from the x-axis and v± 
denotes the normal component of the velocity pointing outwards, then: 



Rw' = -v y cos(0) + Ruj k sin(0) = tan(0) = 

Ruji 



Also 



Ruj 

v± = —v y sin(0) — Rujk cos(0) = — tan(0)i?ojfc sin(0) — Rujk cos(0) 



cos(0) 



implying that 



v' x = v± cos(0) + i?wsin(0) = —Rujk + i?wsin(0), 
v'y = v± sin(0) — Ruj cos(0) = — v y — Ruj cos(0), 

where v' x and v' y are the components of the velocity after the collision with disk Dk+i- 
If we choose v y and Rujk large compared to Ruj, the particle will follow a trajectory 
L' very close to the reversed L on the way from disk D k+1 to disk Dk- By bounding 
the angle of reflection from disk D k +i (choosing v y and Ruj k as large as we need), we 
can ensure that L' hits disk Dk is a small neighborhood of (2k — 1, — R). During the 
flight of the particle to and from the disk Dk+i, we can reset the angular velocity of 
disk Dk to a new value Roj' k with \Roj' k \ large enough so that after the second collision 
with Dk, the particle leaves the system with no additional collisions. By choosing v x , 
v y , Ru>k and Roj' k sufficiently large and large enough compared to w, we can ensure 
that this procedure can be done in an arbitrarily short time. 

Suppose now that uj' ^ 0. Then we have to hit the disk Dk+i at a slightly different 
angular position 0' due to the playground geometry. Then 

„ / ,„ts „ . Ru!kSm(9') — Ruj' 
Ruj' = -v y cos(6') + Ruj k sin(0') => v y = — 



y cos(0') 
v± = —v y sin(0') — Ru>k cos(0') = 

Ruj ksin(9') - Ruj' . , A Rui k - Ruj' sin( 

-( — r ) sin(0 ) — Ruj k cos( 1 



cos(0') > ^ > * v > cos(0') 
Implying that 

v' x = v ± cos(0') + Ruisin(9') = -Ruj k + (Ruj' + Ruj) sin(0') 

v' y = v ± sin(0') - Ruj cos(0') = -v y + (Ruj' - Ruj) cos(0') 

If we choose 0' sufficiently close to 0, v y and Rujk large compared to Ruj and Ruj', we 
can ensure that trajectories from and to disk Dk are very close to the to L. Therefore, on 
the way back, the particle hits disk Dk in a small neighborhood of (2k — 1. —R) and, 
by choosing \Ruj' k \ large enough, we can send the particle out of the system without 
additional collisions. Again, by choosing v x , v y , Ru>k, and \Ru>k \ sufficiently large and 
large enough compared to Ruj and Ruj' , we can do this procedure in an arbitrarily short 
time. □ 
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7. Proof of Proposition [3] 

In this section we complete the proof of Proposition [3] In section [6] we presented a 
construction of a sample path ax from any admissible state X to an arbitrary particle- 
less state Xq € J?o- In subsection 17.11 we would extend this sample to a given state 
Yq e Qq, obtaining a sample path a from X to Yq defined on some time interval 
[0,T]. To finish the proof of Proposition [3] we are also required to show that there 
exists a canonical neighborhood £ of a such that each sample path in £ ends in Uq, a 
neighborhood of Yq. In subsection [772] we treat the situation when the particles in the 
initial state X either collide a disk non-tangentially or exit the system. We finish the 
proof in subsection l7.3l bv treating the remaining situation with tangential collisions. 



7.1. From any state in Qq to any state in Qq. 

Lemma 9. Given Xq, Yq £ [2q and T > 0, there exists a sample path a : Xq — » Yq on 

[o,n 

Proof. Denote the angular positions and velocities of the disks in Xq andYbby (<pi,u)x), 
(ip N ,uj N ) and (tpifU^), • • • , (<p' N ,w' N ) respectively. 
Divide time interval [0, T] into N equal subintervals of length Suppose we 
can set the angular position and velocity of disk Dn to (uj' N ,ip' N ) on time interval 
(0, jf ), where <p' N = (p' N + ^ N ^ T ui' N . Then if we guarantee that no collisions oc- 
cur with disk Dn on time interval [£,T], Djy will have angular position and veloc- 
ity ((p' N ,uj' N ) at time T. Similarly we can proceed with setting angular position and 
velocity of disk -Dj\r-i to (ui' N _ 1 ,ip' N _ 1 ) on time interval (£, where (p' N _ 1 = 
(p' N _ 1 + - — N ^>n-i< ensurm g that if no collisions happen with disk Dat_i on time 
interval [^,T], its angular position and velocity would be (ip' N _ 1 ,ui' N _ 1 ) at time T. 
And so on. In order for this procedure to work, we need the following lemma: 

Lemma 10. Suppose disk Dj rotates with angular velocity uj and there are no particles 
present in the system. Given time t > 0, angular position ip 1 , and angular velocity oj', 
there exists a sequence of particle injections on time interval (0,i) from the left bath 
such that: 

— at time t the disk Dj has the angular position and velocity (tp' ,uj'), 

— at time t all the injected particles have left the system, 

— on time interval (0, t) the injected particles only hit disks D\, • • • , Dj. 

Proof. The result of LemmafTOlis achieved by the application of the Lemma|6]twice: 

Fix some ui > 3+ *" . Apply the Lemma|6]to set the angular velocity of Dj to W\ in 
time |. Let T\ < | be the time of the unique collision with disk Dj. Suppose we wait 
for some time r < | (to be defined later) after | and then apply the Lemma [6] again 
to set the angular velocity of Dj to uj' in time |. Let T2 < | be the time of the unique 
collision with disk Dj counted from the time | + r. 

Then at time t, the angular position and velocity of Dj are 

t 2t 
(Yp + ticj + (- - n +r + t 2 )uj 1 + (— — t — t 2 )uj'] mod l,a/) 
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Let ip = [ip + nu + (| — Ti + T2)cji + (y — T2)<jj'\ mod 1; this is a fixed number 
since all the variables in the expression are fixed. Then we want to pick t < | such that 
p' = (p + t(lui — u))) mod 1. This is not a problem since < | by the choice 

of wi, i.e. in time |, a disk rotating with angular velocity ui — a/ makes full revolution 
and thus, starting at ^, passes through the angular position p' at some time r < |. □ 
This completes the proof of Lemma|9] □ 

7.2. Proof of Proposition^ No Tangential Collisions. In the situation when all parti- 
cles in X either collide with a disk non-tangentially or exit the system under the evo- 
lution of the system with no particle injections, Proposition [3] follows from the lemma 
below: 

Lemma 11. Let a be a sample path from a state X to a state Y on time interval [0, T] 
such that each particle present in the system at any time subinterval of [0, T] follows 
a proper projected particle path. Then for any neighborhood U of Y, there exists a 
canonical neighborhood S of a such that each sample path in £ ends in U. 

Proof of Lemma\TJ\ Denote by c the sequence of injections that generates a. Let £' 
be any canonical neighborhood of er; denote by U' x the neighborhood of X and by C 
the canonical neighborhood of c such that each sample path in S' starts with an initial 
condition in U' x and follows a sequence of injections from C . Define / : U' x X C — > fl 
as follows: if a' G £' is a sample path that starts at state X' G U' x , is generated by 
a sequence of injections c' € C, and ends at state Y', then f(X', e') = Y'. To prove 
lemmafTTl it is enough to show that / is continuous at (X, c). 

We assumed that along sample path a each particle follows a proper projected parti- 
cle path, i.e. it is only allowed to collide with disks non-tangentially. The continuity of 
/ follows from the following facts: 

- If a particle does not collide with dF on time interval [ti,T2], its position and 
velocity change continuously. 

- If a particle collides with the wall on time interval [n, T2] and it is not involved 
in any other collisions with dT on time interval [tj., ra], then its final position and 
velocity depend continuously of its initial position and velocity. [This fact follows 
from the continuity of the billiard flow at collisions] 

- If a particle collides with a disk non-tangentially on [n , tq] and neither the particle 
nor the disk is involved in other collisions on time interval [n, T2], then the parti- 
cle's position and velocity as well as the disk's position and angular velocity depend 
continuously on their initial positions and velocities. [Follows from Lemma|3] 

- If a particle exits through or 7^ on time interval [n, T2] and does not collide 
with d r on time interval [ti , tq\, then the coordinates of the other particles and disks 
are independent from the coordinates of the exiting particle on time interval [ti , 00). 

- The position and velocity of an injected particle depend continuously on the in- 
jected parameters. [Follows from Lemma|2] 

□ 



7.3. Proof of Proposition^ Tangential Collisions. When we constructed a sample path 
ax from X to some particle-less state Xq G J?o in section|6j we first assigned a path in 
r to each particle in X from its initial position to an exit. In order for a particle to follow 
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such a path, the disks had to be set to unique angular velocities at collisions (with Rui = 
v = v t at tangential collisions). Then we showed that setting the disks to nearby angular 
velocities at collisions makes particles follow nearby paths with collisions happening 
at nearby times; this was crucial for choosing appropriate particle paths such that no 
simultaneous collisions with same disks occur. 

When choosing nearby paths in order to avoid simultaneous collisions with same 
disks, we might have to require that at some tangential collisions, R times the angular 
velocities of the disks are not equal to the velocities of the colliding particles. And near 
a tangential collision with Rui ^ v = v t , final position and velocity of a particle does 
not depend continuously on initial position and velocity unlike in the situation with 
Rui = v = Vt- That prevents us from direct extension of Lemma [TT] to the situation 
when tangential collisions might occur. 

However, Proposition|3]only requires us to choose a sample path a and a canonical 
neighborhood E of a such that each sample path in E ends in the given neighborhood 
Uq. By making the size of the discontinuity small enough, we will still be able to ensure 
that every sample path in E ends in Uq. 

The following two Lemmas imply Proposition [3] 

Lemma 12. Given e > 0, there exists time Tx > 0, a sample path ax on [0, Tx] from 
X to some particle-less state Xq and a canonical neighborhood Ex of ax, such that 
each sample path in Ex ends in an e-neighborhood U e (Xo) of Xq. 

Remark 3. Note that state Xq depends on the choice of ax, while the size e of the 
neighborhood U^Xq) around Xq does not. 

Lemma 13. Given Yq £ {2q and a neighborhood Uq ofYQ, there exists e > such that 
for any state Xq G f^o, there exists time Tq > 0, a sample path ao : Xq — > Yo on 
[0,Tq], and a canonical neighborhood Eq of aQ in which each sample path starts in 
t/ e (Xo) and ends in Uq and for any point Y € U £ (Xq), there exists a sample path in 
Eq that starts at Y. 

Lemma [Qlfollows directly from Lemmas l9l and fTTI 

Proof of Lemma \12\ 

As in section|6] denote the initially assigned path in r traced by the j th particle by 
jj. Let t( , • • • , T^yj be the times of collisions with disks, D^y^ , -Dfc(j, n (j')) an d 
let ui{, ■ ■ ■ , ^^j) be the required angular velocities (with Ruj\ = v = vt at tangential 
collisions). 

In the following Lemma, we assume that at time the j th particle is the only particle 
in the system in order to ensure that the system is defined at all times. 

Lemma 14. Assume the j particle has m(j) > tangential collisions before it has 
a non-tangential collision or exits the system. Then there exist open neighborhoods V 3 
of (Qj i v j). If °f ' ' ' , ^n(j) °f u °n(j) suc ^ that for any choice of angular velocities 
(uj{Y € • • • , O^fj))' ^ n,(j)' 'f a P art ' c ^ e starts with a position and a velocity from 
V J , it possibly collides with Dk^j i), ■ ■ ■ ,-Dfc(j.m(j)) set to (uj{)', (uj J m ^)', collides 
with D{ i . /• .\ , i-i , • • • , Di i ■ t y, set to (uj 3 ,■,,,)', ■ ■ ■ , (uP , -s ) , and exits the system. 
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We will prove LemmafPfl after finishing the proof of Lemma IT2l 
Let TL be an upper bound on the time it takes for the j th particle to exit T along 
all possible paths described in Lemma fT4l and let Tx = maxi< j<k{Tx}. Define gj : 
V 3 x I{ x • • • x I° n tj\ — > &o as follows: if a particle starts with a position and a velocity 
from V j , and D k{jA) , ■ ■ ■ ,D h(J ^ {j)) are set to {u)[)' G I{, {u 3 n(J) )' € I 3 n(J) at 
potential collision times, let gj(q, v, (w^)', • • • , (k^ypO be the state of the system at 
time Tx- g 3 is continuous at [^Vq, uj\, ■ ■ ■ ,Lj° n ,^); so there exist sub-neighborhoods 
Vi of V 3 , l{ £ of l{ ■ ■ ■ , f u) e of P n[j) such that for any Vj ) € V 3 , (wj)' G P £ , 



|<fe(<&,«i, K)', ■ ■ • , K w )') - .9,(9°, <4 ■ • • , < e/2. 

Now we are ready to deal with fc-particle system. As in section [6] we can choose 
a path for each particle in X such that upon each disk collision, (uj 3 )' £ I 3 e and no 
simultaneous collisions with same disks occur. That defines ax on [0,Tx]; let X$ be 
the state where ax ends. To define Ex choose further sub-neighborhoods of I 3 e 's to 
ensure that each sample path is Sx is defined up to time Tx - Then each sample path in 
Sx ends in an e-neighborhood U e (Xo) of Xq by the above inequality. □ 

Proof of Lemma \14\ 

By Lemma |7] and the fact the near a tangential collision with Rcu = v — v t , parti- 
cle's final position and velocity depend continuously on its initial position and veloc- 
ity, there exist open neighborhoods Vi{ of v®), I 3 , . N . . „ of u> 3 , ., I 3 , „ 

of u) 3 , ... such that for any choice of angular velocities (uj 3 ,••,,,)' € lL,-\ ■ ■ ■ , 
n(j) •> 6 \ m(j)+l/ ^ m(j)+l' ' 

{uj 3 n ^)' E ^n(j) if a particle starts with a position and a velocity from Vq , it possibly 
collides with D k(hl) ,■• • , A^my)) settowj,--- collides with Z? fc(j:m(j)+1) , 

• • • . D k(j,n(j)) set to ■ ■ • , and exits the system. 

Now we would like to allow an open neighborhood of angular velocities around the 
m th tangential collision. The neighborhood V 3 can be split into two parts: (VJf) c U 
(V^)n C = Vq J , where (Vq) c denotes the set of initial positions and velocities such 
that a particle with a position and a velocity from (Vq) c will collide with -Dj^y m y\\ 

provided that D 3 k{jA) , ■■■ , £>iy, ro y)_i) are set to an S ular velocities wj, • • ■ , ^ n y)_i 
before potential collisions. 

Then the position and velocity of a particle after collision Di-(j, m (j)) depend con- 
tinuously on its initial position and velocity in (Vq) c even if the angular velocity of 
Dk(j,m(jJ) is n °t equal to Also, since particles in (V J ) nc do not collide with 

Dk(j,m(j))7 they exit the system provided the angular velocities of the disks D k ^^, 

• • • , f fc(j - ,m(j)-i) are set to w i' ' ' ' > at appropriate times and angular veloc- 
ities of the disks D k{j ^ m{j)+1) , D k{ , jMl)) are set to values from I 3 m{j)+1Q , 

F n before collisions. 

n(j),0 

Therefore there exists an open neighborhood i^y\ of c*^, .s an d open subneigh- 
borhoods V 3 , of VJ! , I J , , , i -\ of IL., I 3 , * , , of I 3 , ~ such that if 
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a particle starts with a position and a velocity from V^.y it possibly collides with 
D k(j,i)< " ' ' D k(j,(m(j)-i)) with angular velocities uj{, ■ ■ ■ , ^ n(j) _p possibly collides 
with D i(j,m(j)) with an 8 ular velocity {u 3 m(j) )' G I 3 m{j) , collides with D kUMj)+1) , 
D fcO>«) with angular velocities G ^ - )+ i, ro(j) , ■". G 

^"0) anc ' ex ' ts ^ e svstem - The remaining tangential collisions are treated sim- 
ilarly. □ 
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